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Abstract 

It has been shown that a gravitational dual to a superconductor can be obtained by 
coupling anti-de Sitter gravity to a Maxwell field and charged scalar. We review our 
earlier analysis of this theory and extend it in two directions. First, we consider all 
values for the charge of the scalar field. Away from the large charge limit, backreaction 
on the spacetime metric is important. While the qualitative behaviour of the dual 
superconductor is found to be similar for all charges, in the limit of arbitrarily small 
charge a new type of black hole instability is found. We go on to add a perpendicular 
magnetic field B and obtain the London equation and magnetic penetration depth. We 
show that these holographic superconductors are Type II, i.e., starting in a normal 
phase at large B and low temperatures, they develop superconducting droplets as B is 
reduced. 
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1 Introduction 



Holographic superconductors are strongly coupled field theories which undergo a supercon- 
ducting phase transition below a critical temperature, and which have a gravity dual in the 
sense of the AdS/CFT correspondence [b . The reason for considering these systems is that 
they admit a large iV limit in which many aspects of the physics can be studied directly, 
despite the theory being inherently strongly coupled Q The AdS/CFT correspondence ap- 
plied to these systems therefore provides a tractable model for non-standard dynamical 
mechanisms driving superconductivity. 

The existence of holographic superconductors was established in (2J [3] . From the (d 
dimensional) field theory point of view, superconductivity is characterised by the conden- 
sation of a, generically composite, charged operator O for low temperatures T < T c . In the 
dual dimensional) gravitational description of the system, the transition to supercon- 

ductivity is observed as a classical instability of a black hole in anti-de Sitter (AdS) space 
against perturbations by a charged scalar field if). The instability appears when the black 
hole has Hawking temperature T = T c . For lower temperatures the gravitational dual is 
a black hole with a nonvanishing profile for the scalar field tp. The AdS/CFT correspon- 
dence relates the highly quantum dynamics of the 'boundary' operator O to simple classical 
dynamics of the 'bulk' scalar field ip [U |5] . 

In this paper we will expand on previous works by considering two aspects of holographic 
superconductors in d = 3 spacetime dimensions in some depth. Firstly, our previous com- 
putations of the conductivities in the superconducting phase [3] were performed in a limit 
in which the charge of the operator O was taken to be large. This limit had the virtue of 
eliminating the backreation of the scalar field ij) on the spacetime metric. We had expected 
this limit to capture accurately the essential physics and, by considering the theory with 
backreaction onto the metric, we show below that indeed it does. We will furthermore 
discuss the full behaviour of the theory as a function of the charge of O. Remarkably, we 
find that superconductivity persists for arbitrarily small charge. Even when O is neutral, a 
condensate forms at low temperature. This last fact indicates that there are two different 
mechanisms driving the instability, as we will discuss below. 

Backreation on the metric will result in a coupling between the electric and energy 

currents. In addition to electrical conductivity, we also compute thermal and thermoelectric 

1 The large N limit in the best understood cases of AdS/CFT is the large N limit of a gauge theory (with 
order iV 2 ultraviolet degrees of freedom) rather than a vector model (with order N degrees of freedom). 
Unlike vector models, these theories do not become Gaussian in the large N limit. 
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conductivities. Given that our computations are in the clean limit — there are no impurities 
[6] — our system is translationally invariant. This results in the Drude peak becoming a 
divergence in the electrical conductivity at u = 0, even in the normal phase. The divergence 
due to translational invariance (i.e. conservation of momentum) should be distinguished 
from the symmetry breaking infinite superconductivity current, as we discuss in some detail 
below. 

In the second half of the paper we shall study the behaviour of the superconducting 
phase under an external magnetic fieldj^] The theories that we are considering do not have a 
dynamical photon; the U{\) symmetry that is spontaneously broken in the superconducting 
phase is a global symmetry. At zeroth order the theory is perhaps more accurately described 
as a charged superfluid. One might therefore have suspected that little could be said 
about the crucial magnetic properties of superconductors, such as the Meissner effect and 
the difference between type I and type II superconductors. However, we shall show that 
pessimism is unjustified. The key aspect of the Meissner effect involves the generation of 
currents as described by the London equation. We see that this generation indeed occurs 
in holographic superconductors. Furthermore, by weakly gauging the theory we explicitly 
determine the type of the superconductor and compute the mass generated for the photon. 

The symmetry breaking condensates that we study in this paper all have an s wave 
character. It is known that the AdS/CFT correspondence can also describe condensates 
that appear to share many properties of p wave superconductors IT2| [13] . It would be 
interesting to address the range of questions we consider below for the p wave duals. 

Since we are using AdS/CFT to describe superconductivity, it is natural to ask what is 
the connection between superconductivity and conformal field theories. One connection is 
that the CFT describes a quantum critical point, i.e., a phase transition at zero temperature. 
In some cases the large fluctuations associated with the quantum phase transition can induce 
the pairing responsible for superconductivity [Tl] ■ The quantum critical theory is the theory 
of these fluctuations, and AdS/CFT will provide us with models for such theories. In this 
work however, we simply use the fact that the AdS/CFT correspondence has been extended 
to nonconformal field theories as well. We will break the conformal invariance by adding a 
background charge density. 

This paper is organized as follows. In the next section we write down the equations 

2 There have been a couple of previous discussions of holographic superconductors in the presence of 
magnetic fields [7J [5]. These have not settled the question of the type of superconductor. Some early 
attempts to add magnetic fields are [§1 IIP]. 
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for the gravity dual of a superconductor. The rest of the paper is devoted to investigating 
the solutions to these equations. In section three we study the static solutions describing 
hairy black holes which are dual to the superconducting phase with nonzero condensate. 
The following section contains a discussion of transport phenomena in the superconductor, 
using perturbations of the black hole. In section five, we introduce an orthogonal magnetic 
field and show that the superconductor is Type II. The next section contains a discussion 
of the currents induced by the magnetic field, and shows how the London equation is 
recovered. In section seven we explicitly weakly gauge the theory and compute the photon 
mass in the superconducting phase. Some comments on the relation and difference between 
the holographic approach to superconductivity and Landau-Ginzburg theory are in section 
eight. We conclude with a short summary and some open questions. 

2 The bulk equations for a holographic superconductor 

The best understood examples of the AdS/CFT correspondence involve AdS spaces that 
are part of a full ten or eleven dimensional solution to string or M theory. The low energy 
fields that propagate in the AdS space are obtained as consistent truncations of the higher 
dimensional theory. However, in the field of 'applied AdS/CFT' (one is thinking usually 
of applications to QCD or condensed matter physics) a different philosophy is possible. A 
phenomenological approach is taken in which the classical fields which propagate in the bulk 
and their interactions are chosen by hand to capture the physics of interest. The two main 
drawbacks of this approach are that firstly it is not possible to compute many quantities 
away from the large N limit (because the bulk theory is likely not part of a consistent 
theory of quantum gravity) and secondly we do not have an explicit description of the 
field content and interactions of the dual field theory. We believe it will be possible and 
certainly interesting to realise holographic superconductors as truncations of string theory. 
(A concrete suggestion for embedding p wave superconductors was made in [13].) For the 
moment we will work with a phenomenological model. 

What are the minimal ingredients we need to describe a holographic superconductor? 
We are interested in the continuum limit, as AdS/CFT has not yet been developed for 
lattices, and therefore our field theory will have a conserved energy momentum tensor T^ v . 
For simplicity we will be working with a theory that is Lorentz invariant at high energies 
and so the indices /x, v run over t,x,y. The AdS/CFT correspondence has recently been 
generalised to non-relativistic theories |15[ I16j . and it would certainly be interesting to 
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adapt our study to those cases. Furthermore we need a global U{\) symmetry in the field 
theory, and therefore we will have a conserved current J^. Finally, as we wish to break this 
U{\) symmetry spontaneously, we need a charged operator O, which will condense at low 
temperature. 

The most basic entries in the AdS / CFT dictionary jH [5] tell us that there is a mapping 
between field theory operators and fields in the bulk. In particular, T^ v will be dual to the 
bulk metric g ab , the current will be dual to a Maxwell field in the bulk, A a , whereas the 
operator O will be dual to a charged scalar field ip (which is therefore necessarily complex) . 
Here a, b run over the four bulk coordinates t, x, y, r. The next step is to write down a 
minimal Lagrangian involving these fields. The Lagrangian density for a Maxwell field and 
a charged complex scalar field coupled to gravity is 

C = R + A _ ±F ab F ah - V{\^\) - \Vi; - iqA^ 2 . (2.1) 

As usual we are writing F = dA. Most of our work will revolve around solving the equations 
of motion that are obtained from this Lagrangian. These are the scalar equation 

- (V„ - iqA a ) (V a - iqA a ) + ~V'M) = , (2.2) 

2 

Maxwell's equations 

V a F ab = iq [4>*(V b - iqA h )^ - ^{V b + iqA h )^*] , (2.3) 

and Einstein's equations 

Rgb _9a*R_3to = l_ FacFb c _ 9_^ F c dFcd _ 9_^ vm 

-M| V ^ - iqAiP\ 2 + i [(V ^ - iqA a iP)(V b r + iqA b r) + a ~ b] . (2.4) 

Note in these equations that q is the charge of the scalar field. 

The first step will be to find static black hole solutions to these equations. These 
solutions describe the equilibrium phases of the theory. It is in these solutions that we see 
a phase transition as the temperature is lowered. Once we have the equilibrium solutions 
we can look at perturbations away from equilibrium, which describe transport processes. 



3 Normal and superconducting phases 
3.1 The hairy black hole ansatz 

In the AdS /CFT correspondence, plasma or fluid-like phases of the field theory at nonzero 
temperature are described by black hole solutions to the bulk gravitational action [T7] . 



5 



To study conductivity and other transport properties of a charged plasma, the first step 



therefore is to find black hole solutions to our theory (2.1). Superfluidity and superconduc- 
tivity are associated with symmetry breaking, and thus we search for solutions in which the 
charged scalar field has a nontrivial expectation value. 

One further ingredient is necessary. The simplest theories described by the AdS/CFT 
correspondence are conformally (in particular, scale) invariant. In a conformal field theory 
in Minkowski space, such as ours, in the absence of another scale, all nonzero temperatures 
are equivalent. Therefore, if we wish to obtain phase transitions at a critical temperature, 
we need to introduce another scale. Among the various scales to choose from, we introduce 
a finite charge density p (equivalently, a finite chemical potential p). In 2+1 spacetime 
dimensions, p has dimensions of mass squared whereas p has mass dimension one. Our 
primary motivation for introducing a scale in this manner is the observation in j2] that it 
leads to a superconducting instability at low temperatures. If we were to compare with 
real experimental systems, our model describes a quantum critical theory that has been 
deformed by doping [18 . Another possible connection is to relativistic systems such as 
graphene held at finite gate voltage. 

As we shall review below, a charge density p in the system corresponds to giving an 
electric charge to the black hole. We shall consider magnetic charges also in a later section. 
The upshot is that we are looking for electrically charged plane-symmetric hairy black hole 
solutions. Thus we take the metric ansatz 

ds 2 = -g(r)e-^dt 2 + ^ + r 2 (dx 2 + dy 2 ) , (3.1) 

g(r) 



together with 



<t>{r)dt, 7p = ^( r ). (3.2) 



We now look for solutions to the above equations of motion with this form. 

It is immediately seen that the r component of Maxwell's equations implies that the 
phase of ip must be constant. Without loss of generality we therefore take ifi to be real for 
the background. The scalar equation becomes 



+ ( *- - % + 2 ) * 4- - 1 V W = , (3.3) 

g 2 r ) g l 2g 



Maxwell's equations become 
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while the tt and rr components of Einstein's equations yield 

x < + rr + "££pl =0 , ( 3. 5) 

9 

H — * 1 \~ ~9 F9 H — o I — = - I 3 - 6 ) 

2 4g gr r 2 gL 2 2g 2g 2 

Note that both these equations only have first derivatives, although they appear squared. 
The xx component of Einstein's equations is not independent and follows from differenti- 
ating the two equations above. We will specialise for concreteness to the simple potential 

VM = ~1?. (3.7) 

This is the conformal mass term for a scalar in AdS^ and is above the Breitenlohner- 
Freedman bound for stability. We choose this value following [3 j . The effect of other masses 
is discussed in [19] . A stringy embedding of this model would of course fix the full potential. 

If one takes the limit q — > oo keeping qip and q<f> fixed, the matter sources drop out 
of Einstein's equations (3.5) and (3.6), while the scalar and Maxwell equations ( 3.3|3.4 ) 
remain essentially unchanged. This is the probe limit studied in [3 . Our first objective in 
this paper is to go beyond the probe limit. We will solve the full set of equations (with 
finite q) numerically by integrating out from the horizon to infinity. By considering a series 
solution at the horizon — the horizon radius r + is defined through the requirement that 
g(r+) = — one finds that there are four independent parameters at the horizon 

r+ , V+ = i/>( r +) > E + = c t>'(r + ) , x+ = X(r+) ■ (3.8) 

The third of these quantities is the value of the electric field at the horizon. The scalar 
potential <p itself must go to zero at the horizon in order for the gauge connection to 
be regular. These quantities determine the Hawking temperature of the black hole (for 
instance, from regularity of the Euclidean solution) 

T = ((12 + 4^) e -X+/2 _ J2 E 2 + eX+ /2^j _[±_ _ (3 9) 

At infinity we have the following parameters that determine the charges of the black hole 
and the expectation values of scalar fields. The charge density, p, and chemical potential, 
fi, are read off [H [5] (more explicitly, see for instance [20J) from the asymptotic value of the 
scalar potential as r — > oo 

= /i- £ + •••. (3.10) 
r 

The general asymptotic behaviour of the scalar field as r — ► oo is 

1>= + ^2- + --- ■ (3-11) 



This simple falloffis another reason for choosing the mass (3.7). This value of the mass falls 
within the range in which there is a choice of admissible boundary conditions at large radius 
|21j . Depending on the choice of boundary conditions, we can read off the expectation value 
of an operator O2, of mass dimension two, or of an operator Oi, of mass dimension one. 
Specifically, for a stable theory we must either impose 

= , and (0 2 ) = v / 2V (2) , (3.12) 

or 

V> (2) = , and (d) = V2xfj^ . (3.13) 

The factor of y/2 is following [3] and is a convenient normalisation. 

In order for the Hawking temperature of the black hole to be the temperature of the 
boundary field theory, we must impose 

X^O, as r — ► 00 . (3-14) 

This is a statement about the normalisation of the time coordinate t relative to the grav- 
itational redshift, as determined by the normalisation of the r coordinate. In practice we 
can implement this boundary condition by taking an arbitrary x+ at the horizon, obtaining 
the asymptotic value of Xi an d then rescaling time t — > at, in order to set the asymptotic 
value to zero. Said another way, we are using the scaling symmetry of the metric, gauge 
field, and equations of motion 

e x — > a 2 e x , t — > at , (j)— >■(/>/ a, (3.15) 

to set x = at the boundary. 

There are furthermore two scaling symmetries of the equations of motion that we can 
use to set L = 1 and r+ = 1 when performing numerics. The first is 

r — ► ar , t — ► at , L — ► aL , q^q/a, (3.16) 

which rescales the metric by a 2 and A = <fi dt by a. The second is 

r^ar, (t,x,y) (t,x,y)/a , g^a 2 g, 4> , (3.17) 

which leaves the metric and A unchanged. After these scaling actions we are left with one 
parameter in the Lagrangian which is physical, the charge of the scalar field q, and two 
parameters that determine the initial data at the horizon, ip + and E + . Integrating out 
from the horizon to infinity gives a map 

(^ + ,ii; + )^(^,p,V (1) ^ (2) ,e). (3.18) 



8 



We have included the mass of the black hole solution on the right hand side of this expression. 
We have denoted the mass by e, as it is to be interpreted as the energy density of the field 
theory. The mass is to be read off from the large r behaviour of g and X- Assuming that 
X —* as r — > oo we have 

,2 ^(1))2 _ eL 2 /2 + 4^,(1)^,(2) /3L 2 
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r 
1? 

V ^ 

L 2 " 



+ 



2L 2 
eL 2 



+ 



+ 



2r 



+ 



(3.19) 
(3.20) 



Finally, upon imposing either of the boundary conditions (3.12) or (3.13), the map (3.18) 
reduces to a one parameter family of solutions for each value of the scalar field charge q. 
We can think of this parameter as being the temperature of the theory at a fixed charge 
density. Depending on whether the scalar ip is nonzero or not for this solution, we will be 
in the superconducting or normal phases, respectively. 



3.2 Phase diagram 

In figure [T] we show the condensate for the charged scalar field as a function of temperature 
with the charge density held fixed. We plot q{0) since this is the quantity which is finite 
in the probe (large q) limit. This plot has been obtained by numerically solving the differ- 
ential equations in the previous subsection using a shooting method. Although there is an 
ambiguity at this point in the normalisation of the scalar field, we will shortly relate the 
condensate, in our normalisation, to physical quantities such as the gap in the frequency 
dependent conductivity. 

The most important feature of the plots is that in all cases there is a critical temperature 
T c below which a charged condensate forms. This is the symmetry breaking phase transition 
to a superconducting phase. For T > T c the solution is simply the Reissner-Nordstrom-AdS 
black hole. That is (with L = 1) 

*=*=»• »- 2 -K ri+ £) + £' (3 - 2i) 

At these high temperatures, there are no hairy back hole solutions. At T = T c the Reissner- 
Nordstrom-AdS solution becomes unstable against perturbations of the scalar field. As 
pointed out in [2] , this instability can be understood directly from the fact that the coupling 
of the scalar to the gauge field through covariant derivatives induces an effective negative 
mass term for the scalar field. This term becomes more important as the temperature is 
lowered at fixed charge density, eventually driving the scalar field tachyonic. For T < T c 
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we find that hairy black hole solutions do exist and have a lower free energy than the 'bald' 
black hole. 

The general form of the curves is similar to the well known case of, for instance, BCS 
theory. The condensate turns on at T = T c following a square root law (Oi) ~ T*(l — 
T/Tc) 1 / 2 , as is typical for mean field theory treatments of second order transitions. As the 
temperature is taken to zero, the condensate tends to a finite value in terms of the scale 
set by T c . It is interesting to note the dependence of the condensate on the charge of the 
operator O. The most striking effect is seen in the theory with operator 0\, in figure [T^,. It 
had been found in |3j that in the probe limit (q — * oo) the condensate appeared to diverge 
as T/T c — ► 0, perhaps indicating an instability of the theory. In figurejl^ we can see how the 
full backreacting system cures this divergence. At smaller values of the charge q there is no 
sign of a divergence. As q is increased the plot starts to curve upwards as we approach low 
temperatures and appears similar to that of the probe limit. Unfortunately, it is difficult to 
get the numerics reliably down to very low temperatures so we have not been able to see the 
resolution of the divergence at large but finite q explicitly. Clearly it would be extremely 
desirable to have a more direct approach to the zero temperature properties of holographic 
superconductors. One numerical observation we can make is that the electric field on the 
horizon, E+, appears to go to zero at low temperatures. Note that, in contrast, for the O2 
theory the condensate always stays bounded as T/T c — ► and probe limit is approached 
from above. We believe that this behavior is more typical of holographic superconductors. 

As was noted in [3] a second order transition is only possible in our 2+1 dimensional 
system at finite temperature because we are working in a large N limit. The large N limit 
suppresses fluctuations of the fields, in particular, the massless fluctuations associated to 
the spontaneous breaking of our global £7(1) symmetry. Away from the large N limit, these 
fluctuations will lead to infrared divergences which will destroy the long range order (the 
Coleman-Mermin- Wagner theorem) in the low temperature phase. It would be extremely 
interesting to capture these fluctuations within an AdS/CFT framework and perhaps to 
exhibit an algebraic order with spatially separated correlations falling off like (Ax) _1//Ar# . 
A discussion of large iV expansions and symmetry breaking in two dimensions can be found, 
for instance, in [22]. A related question is whether any analogue of the Berezinski-Kosterlitz- 
Thouless transition can be seen in this or other AdS/CFT systems. Finally, in section 
seven below we will gauge the U(l) symmetry (with 2+1 dimensional photons) and see that 
the Goldstone boson is eaten by the photon, which becomes massive. Once there are no 
Goldstone bosons, there are no longer IR divergences. 
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a) T c b) Tc 

Figure 1: The value of the condensate as a function of temperature, with the charge density 
held fixed, for the two different boundary conditions: a) from bottom to top, the various 
curves correspond to q = 1, 3, 6, and 12; b) from top to bottom, the curves correspond to 
q = 3, 6, and 12. The value q = 1 gives a much larger condensate in this case, achieving 
\J q(02) /T c ~ 21 so we have not plotted it. Note that the large q limit is approached in 
opposite directions in the two cases. 

The critical temperature T c appearing in figure [T] is set by the only other dimensionful 
scale in the system, the charge density p. (In the grand canonical ensemble, we could 
alternatively consider the scale to be set by the chemical potential p,.) Dimensional analysis 
implies that we will have T c oc ^fp. However, the constant of proportionality will also 
depend on the charge q of the operator O. This dependence is shown in figure 2. 

In figure [2] we have also included the value of T c in the probe (large q) limit. This was 
computed in [3] where it was found that T c oc y/qp. The effects of backreaction produce two 
changes with respect to the probe limit. Firstly, T c is suppressed compared to the probe 
estimate everywhere except at very small q. This suppression relative to the probe limit 
can be understood in part from the correction to the Hawking temperature of a black hole 
due to the electric charge: T ~ 12r^ — p 2 . Away from the probe limit, the electric charge 
back reacts on the metric, and this charge decreases slightly the temperature. Secondly, 
T c remains nonzero for all q, including q = oj^] In other words, even neutral operators can 
condense at low temperature. Lest one doubt our numerics, we give a proof of this fact in 
the appendix for the case (0i) 7^ 0. 

3 If T c were to go to zero at some finite q, we would have a zero temperature quantum phase transition 
at that point. We expect such a quantum critical point for more positive masses. 
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'0 1 2 3 4 02468 10 12 

1 b) 9 

Figure 2: The blue solid line is the critical temperature as a function of q. The dashed 
red line is the probe limit (naively extrapolated to all q). a) The dimension one case where 
the probe limit corresponds to 0.2255^^; b) the dimension two case with the probe limit 
0.1173^/g. 

3.3 A new type of instability 

At first sight it is very surprising that a near extremal charged black hole is unstable to 
forming neutral scalar hair. As mentioned earlier, the reason charged scalar hair is expected 
is that the coupling of the scalar to the gauge field through covariant derivatives induces an 
effective negative mass term for the scalar field. This term becomes more important as the 
temperature is lowered at fixed charge density, eventually driving the scalar field tachyonic. 
This mechanism does not apply to neutral scalar fields, so the origin of the hair in this case 
must be qualitatively different. 

The best explanation seems to be the following]^] First we note that although T c is not 
zero at q = 0, it is small. Therefore we will think of these unstable black holes as being near 
extremal. An extremal Reissner-Nordstrom AdS black hole has a near horizon geometry 
which is v4(iS2xl 2 . Our neutral scalar has m 2 = —2/L 2 . The Breitenlohner-Freedman (BF) 
bound governing stability of scalar fields in AdSa is m 2 BF = —(D — 1) 2 /4L 2 . So while our 
scalar field is above the BF bound for AdS 4, it is below the BF bound for AdS%- The AdS2 
radius of curvature is actually smaller than AdS& (L?, = Lf/6 for an extremal black hole), 
but even taking this into account, our scalar is below the BF bound in the near horizon 
region of the black hole. This argument suggests that a Reissner-Nordstrom AdS black hole, 
when coupled to a neutral scalar with m 2 = —2/L 2 , becomes unstable near extremality. A 
proof of this statement is given in Appendix [A} The instability produces the hairy black 
4 We thank M. Roberts for suggesting this. 
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holes we see numerically^] (This argument alone does not explain why the dimension two 
case with faster fall off should be more stable, which is what is seen numerically.) Note that 
extremal Reissner-Nordstrom AdS black holes are not supersymmetric, so their instability 
does not lead to a contradiction with general stability results. 

The physical point to take away from these observations is that there are (at least) 
two distinct physical mechanisms leading to superconductivity in our system. At very 
large charges for the scalar fields, it is the (bulk) gauge covariant derivatives that enhance 
the effective negative mass. At very small charges, it is the fact that the near extremal 
charge of the black hole produces a throat in which even a neutral scalar with sufficiently 
negative mass squared becomes unstable. The crossover between these effects presumably 
corresponds to the crossover visible in figure 2. Although the crossover appears to be 
smooth, one might bear in mind the possibility of a phase transition as the type of instability 
swaps. It would be interesting, of course, to re-interpret this crossover from the dual field 
theory. 

3.4 The hairy black hole action 

We conclude this section by computing the Euclidean action for our hairy black hole: 

S E = - f d 4 x^g^C , (3.22) 



where C is given in (2.1 ) and gg is the determinant of the bulk metric. We first show that, 



when evaluated on a solution, this action reduces to a simple surface term at infinity. From 



the symmetries of the solution (3.1 3.2 ), the xx component of the stress energy tensor only 



has a contribution from the terms proportional to the metric. Thus, Einstein's equation 



(2.4) implies that the Einstein tensor satisfies 

1 



This implies 



or 



G xx = ^r\C - R) . (3.23) 



R = G a a = G t t + G r r + C-R, (3.24) 



C = -G\ - G\ = -1 [(rg)' + foe"*)'^] . (3.25) 



The Euclidean action is then a total derivative 



S E = J d 3 x ^ °° dr[2rge~ x/2 }' . (3.26) 



An earlier example of a charged AdS black hole with neutral scalar hair was given in [23] . That example 
also involved a scalar with m 2 = —2/L 2 ; however it required a black hole with a hyperbolic horizon. 
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The surface term on the horizon vanishes since g(r+) = 0. So we get just the surface term 
at 



S E 



d x 2rge 



-X/2 



(3.27) 



This action diverges as — » oo and must be regulated. The counter terms we need to 
regulate it are standard (see for example [21] )• We require a Gibbons-Hawking term and a 
boundary cosmological constant: 



Set. = f d 3 x^g^(-2K + 4/L) 



(3.28) 



where g^ is the induced metric on the boundary r = and K = g^S/^n u is the trace of 
the extrinsic curvature (n M is the outward pointing unit normal vector to the boundary). We 
also require a term quadratic in the scalar field that depends on which boundary condition 
we choose for ip. If we fix the value of on the boundary, we must add 



Si 



d 3 x^-g OQ ip 2 /L 



(3.29) 



while if we fix the value of on the boundary, we need also an analog of a Gibbons- 
Hawking term, 



So 



d 3 xy/=g^ (2V n»du.il> + ^ 2 /L) 



(3.30) 



The sum Se = Se + S c _t, + Si,2 is now finite in the limit — > cxd. The regularised 
action becomes, combining the cases of the two possible boundary conditions above 



TS E 



d 2 x 



+ 7V> (1) V> {2) 



EL 2 



+ 7V2^ (1 V (2) , 



(3.31) 



with 7 = 2/3 in the first case and 7 = —4/3 in the second case. Here we used the fact that 
the length of the Euclidean thermal circle is 1/T and in the last equality we assumed that 
the ip® were constant in space. For our solutions, at least one of the ip® = 0. We will 
see in Section [5] below that the result (3.31) for the background (equilibrium) configuration 
matches the general expectation for a 2+1 CFT that the grand canonical potential function 
is n = -E/2. 



4 Conductivities 

In this section we will study transport phenomena in our holographic superconductors]^] In 
particular, we obtain the electric, thermal and thermoelectric conductivities as a function 



We set L = 1 for simplicity in the remainder of this paper. This choice only affects the overall normali- 
sation of the bulk action. 
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of frequency. Transport describes the response of the system to small external sources. 
Therefore, we will need to compute the retarded Greens functions for the electric and heat 
currents. In the AdS/CFT correspondence, these correlation functions are computed by 
looking at the linear response of the system to fluctuations of the fields A x and g± x in 
the bulk. These fluctuations are dual to the electric current J x and energy current T tx 
operators in the CFT. At zero spatial momentum, these two fluctuations do not source 
any other modes of the metric, Maxwell field or scalar field. This decoupling simplifies 
considerably the computation, so we shall restrict to zero spatial momentum in this paper. 



4.1 Formulae for the conductivities 

Assuming a time dependence of the form e~ lujt , we linearise the Maxwell and Einstein 
equations above to yield equations governing perturbations of A x and gt x . We start with 
Maxwell's equation 



A" + 



X_ 
2 



A' + 



2q 2 ^ 2 



A a 



9 



-9tx + ~9tx 
r 



(4.1) 



.9" 9 

We can caricature this equation in a way that makes the physics of the conductivity much 
clearer. If we assume that the radial dependence is not essential, then this equation describes 
a photon of mass proportional to q 2 t/j 2 coupled to a metric fluctuation in gt x . Furthermore, 
in the probe limit, we assume the gauge field does not back react on the metric, and we set 
9tx = 0. In this limit, we have at least morally the Higgs mechanism and thus expect the 
charge current response to a magnetic field and the infinite DC conductivity typical of a 
superconductor. If we move away from the probe limit and allow gt x to fluctuate, we find 
remarkably that gt x is governed by only a first order differential equation in the case of zero 
spatial momentum: 



9tx 



-9tx + 4> A x = . 



We can substitute this Einstein equation into (4.1) to find 

,2 



X_ 
2 



A' + 



fj2 



2q 2 ifj 2 



A x = 0. 



(4.2) 



(4.3) 



.9 9 J 9 

The effect of allowing the metric to fluctuate has been to add another contribution to the 

effective mass that scales as 4>' 2 . This extra contribution is the result of restoring translation 

invariance and will lead to the infinite DC conductivity typical of translationally invariant 

charged mediaj^] However, we do not expect this extra mass to affect the response of the 

7 From a quasiparticle point of view, the charged particles of the system are accelerated by the external 
field. There may be scattering events, but translation invariance and a net charge means the end result will 
be acceleration of the entire system. 
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system to magnetic fields. In support of this second claim, we find that if we look at 
fluctuations with nonzero spatial momentum, in addition to having to consider many more 
modes of the metric and gauge field, g% x will satisfy a second order differential equation. 
Therefore we would not be able to simply eliminate gtx from the equation of motion for A x . 



These fluctuation equations (4.2) and ( |4.3[ ) are solved at the linearised level by using 
the equations for the background in section [3j The asymptotic large r behaviour of the 
perturbations is 

A x = A^ + ^- + ..., <7 i , = r 2 ^ + ^ + .... (4.4) 

As is standard in AdS/CFT, the leading term determines a source in the dual theory, while 
the 'normalisable' term will give the expectation value of the dual current. We shall see 
this explicitly shortly. 

Preparatory to calculating two-point correlation functions of the currents, we evaluate 
the quadratic action for perturbations about a solution to the equations of motion. The on 
shell (Lorentzian) action is 

So.s. = dr Jd 3 x ^/=^£ . (4.5) 
Expanding the action to quadratic order in A x and gt x , and using the equations of motion 



(4.2) and (4.3), we find that the quadratic action reduces to a surface term: 

= J d 3 x e*/ 2 (- 9 -e-*A x A! x - g tx9 ' tx + \ {t - X ') g^ | . (4.6) 

There is no contribution from the lower limit of the radial integral since g and gt x both 
vanish at the horizon r = r + . We would like eventually to take the upper limit of the 
radial integral to correspond to the conformal boundary of our asymptotically AdS space, 

— > oo. As usual, S . s . is not finite in this limit and must be regulated by counter terms. 

The counter terms we need to regulate So.s. are the same as those used in section (3.4). 
Adding the appropriate terms together we find that, evaluated on a solution, the regularised 
quadratic action becomes 



s[ 2) = lim (s® + S c . t . + Si 

Too— »00 \ , 

= j d * x (i4 0) 4 1} - 3<£%2 } - (| - ^ (1 ¥ 2) ) AL 0) ) , (4.7) 



and 



Si 2) = lim ( Si 2 ) + Set. + S 2 



roo— >00 



/ d* X - ZgfJgfJ - (| + j^V 2) ) <&2?) • (4-8) 



16 



These regulated expressions are now manifestly finite in the limit Tqq — > oo. 

Equipped with the quadratic action we can obtain the conductivities as follows. Firstly 
note that the mixing of the Maxwell perturbation with the metric mode means that we 
must consider thermal and electric transport jointly. These phenomena are described by 
the matrix of conductivities 





E x 



-{V X T)/T 



(4.9) 



Here J x is the electric current and Q x = Tt x — liJ x is the heat current. We will have all 
currents moving in the x direction and all sources pointing in that direction. The electrical 
conductivity is a, the thermoelectric conductivity is a and the thermal conductivity is R. 
The external fields are an electric field E x and a thermal gradient V X T. Strictly, the currents 
are expectation values. We will drop the angled brackets for notational convenience. The 



matrix appearing in (4.9) is symmetric due to time reversal invariance. 



At a nonzero momentum there will also be mixing with the condensate, which will result 



in more terms in (4.9). We shall not consider these effects here. 



We can solve the equation for the metric perturbation (4.2) to obtain 



(0) , f°° $A X 



9t* = r 2 ( 9 r x >+ I ^dr). (1.10) 
This shows that g^} is independent of the Maxwell perturbation, whereas g^) is completely 



determined by A x ° ■ In fact, by expanding (|4.10|) at large radius and using (|3.10|) and (|4.4l) 



we have 



*£? = 1 4 0) • (4-11) 



Given a solution to the equations, we compute the currents by differentiating the action 
with respect to the boundary values of the dual bulk fields 

5A X 5gl x > 

In general we can study these currents in the presence of arbitrary sources A x ^ and g^) . 
Firstly let us set g^ = 0. Thus there is no source for heat flow, so that V ' X T = 0. In this 
case, the electric and heat currents are given by 

J x = 4 1 ) , Q x = -3^ - = -P4 0) " m4 X) • (4-13) 

In performing the differentiation, we can note that because we are looking at linearised 
equations, A x ^ will be proportional to A x °\ 
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Setting V X T = in (4.9) we can therefore read off the electrical conductivity 



and the thermoelectric conductivity 

Ta 



Err. 



Qx 

Err 



%A X ^ 



■ip 



(4.14) 



(4.15) 



In these expressions we used the fact that Ax is the boundary background potential, and 
that E x - 



■dfAx ■ The relation (4.15) between thermoelectric and electrical conductivities 
is the same as that found in [25], in the absence of a charged condensate. The simple 
relationship between electric and thermoelectric conductivities means that we can focus on 
computing the electric conductivity in the remainder of the paper. 

To obtain the thermal conductivity, we now set A x ^ = 0. Using the fact that our 



backgrounds have either = or t// 2 ) = 0, it follows from (4.12) that 

Qx 



(0) 



(4.16) 



Now we need to relate gi x to a thermal gradient. This is a straightforward computation. 
(See for instance the appendix of [6].) One obtains 



o (0) 
9tx 



iujT 



Combining the last two equations and (4.9) with E x = gives 



le 
u/T 



(4.17) 



(4.18) 



The divergence as uj — > is exactly what we should expect from a translationally invariant 
system. Conservation of momentum means that a DC energy current, which is also a 
momentum, cannot relax. 



4.2 Numerical results for the conductivity 

We make several observations about figures [3] and [4] which display our numerical results for 
the conductivity as a function of frequency to. To understand these plots, it is useful to put 
them in the context of conductivity results from previous work [26] , [25] , [6] , and |3j . 

From the earliest of these papers |26j . we know that the conductivity at vanishing 
charge density and vanishing condensate is a constant independent of u. This independence 
can be understood as a consequence of classical electromagnetic self-duality of the dual 
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Re(cr) 




Re(cr) 
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Figure 3: The dashed red line is the real part of the conductivity at T = T c (for q = 3). The 
blue lines are the same conductivities at successively lower temperature: a) The dimension 
one operator with T/T c = 0.810, 0.455 and 0.201; b) the dimension two operator with 
T/T c = 0.651 and 0.304. There is a delta function at the origin in all cases. 



gravitational theory at the quadratic level. From the next paper [25] in this sequence, 
we have results for the conductivity as a function of charge density in the absence of a 
scalar condensate. If we work in the limit where the charge density is small compared to 
the temperature, we recover the frequency independent result of [26], but in general the 
dependence on lo is more complicated. In particular, we see the minimum in Re(cr) at 
lo = displayed by the dashed curves in figure [3| The Im(cr), not plotted, has a pole at 
lo = 0. From the Kramers-Kronig relations, which follow from causality, one concludes that 
the real part of the conductivity has a Dirac delta function at lo = that is invisible to the 
numerics because of its infinitesimal width. Recall that one of the Kramers-Kronig relations 
is 



Im[cj(a;)] 



1 



-V 



Re[o-(u/)]cL/ 



(4.19) 



vr j_r^ lo' — lo 

From this formula we can see that the real part of the conductivity contains a delta function 
Re[<7(a;)] = tt5(lo), if and only if the imaginary part has a pole, Im[<r(a;)] = 1/lo. There is 
also a Ferrell-Glover-Tinkham sum rule which follows from similar arguments and which 
states that J Re(o")du> is a constant independent of the temperature. Thus the dip in the 
real part of the conductivity at lo = is related to the residue of the pole in the Im(a") and 
strength of the Dirac delta function in the Re(cr). 

This Dirac delta function for T > T c is naively surprising because it implies an infinite 
DC conductivity in the normal phase. This infinite conductivity is not superconductivity 
and results instead from translation invariance. A translationally invariant, charged system 
does not have a finite DC conductivity because application of an electric field will cause 
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0.0 0.5 1.0 1.5 2.0 

w 

^q(0 2 ) 

Figure 4: We plot the real part of the conductivity as a function frequency normalized 
by the condensate, either q{Oi) or \J q(02) as appropriate. This data was taken at low 
temperature, T = 0.03 q(0\) and T = 0.03a/ ^(02) for a variety of charges q = 1,3,6 and 
12. The curves with steeper slope correspond to larger q. There is a delta function at the 
origin. 

uniform acceleration. If we were to break this translation invariance by for example intro- 
ducing impurities, the delta function at uj = would acquire a width for T > T c and the 
conductivity would become finite. This effect of impurities was studied in an AdS/CFT 
setting in [6] . 

In our previous paper [3], we did not see this infinite conductivity above T > T c . 
The reason we did not see it is that we worked in a probe limit where the gravitational 
background was fixed and the abelian-Higgs sectoij^] decoupled. By fixing the background, 
we implicitly broke translation invariance. Technically this occurs because the electric and 
energy currents decouple, as we discussed at the start of this section. For T > T c we 
had a pure Schwarzschild-AdS background and thus recovered the frequency independent 
conductivity of [26 . It was only for T <T C that the Im(cr) developed a pole. 

With this review of previous results, the structure of figure [3] should be clear. We see 

from the dashed curves in figure [3] that Re(cr) has a minimum at uj = for T = T c and 

indeed, although not plotted, also for T > T c since the background becomes the electrically 

charged black hole studied in [25 . If we were to plot Im(cr), we would see a pole in u> = 0, 

and by the Kramers-Kronig relations would conclude that there is also a Dirac delta function 

in Re(cr). In other words, we have infinite DC conductivity for temperatures above T c . 

8 Our bulk matter is like a traditional abelian-Higgs model, but it does not have the usual ip 4 term in the 
potential. 
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For T < T c , this minimum in the Re(cr) at to = becomes increasingly pronounced 
and eventually develops into a gap, as was seen in [3] in the probe limit. For lower q 
the low temperature gap becomes less pronounced, see figure |4j As far as the accuracy of 
our numerics permits, the conductivity still appears to vanish at zero temperature over a 
finite range of small frequencies. The residue of the pole in Im(<r) and the strength of the 
Dirac delta function become much larger as well. There is an additional contribution to the 
strength of the Dirac delta function coming from condensation of the scalar. Note that the 
low temperature plot of figure [4] gives strong evidence that for relatively large values of the 
charge of the scalar field q > 3, the size of the zero temperature gap, uj g , can be associated 
in the dimension one case with q(0\) and in the dimension two case with a/ q(02)- Let us 
take this observation as the definition of tu g . Thus the vertical axis in our figure [T] can be 
thought of as ui g /T c . 

The pole in the conductivity for T < T c is (partly) due to different physics than that 
for T > T c . It is not exclusively due to translational invariance and would persist in the 
presence of impurities. It is instead closely tied to the spontaneous breaking of the U(l) 
symmetry and the expulsion of magnetic fields. In sections [6] and [7] below, we shall discuss 
the physics of magnetic fields and Higgsing in this model in greater detail. 

As we lower T past T c we therefore expect to see a non-analytic change in the strength 

of the delta function, reflecting the onset of superconductivity. Indeed, there is a jump in 

the derivative of the strength of the delta function with respect to temperature, indicative 

of a second order phase transition, as we now discuss. Above the transition temperature, 

we know from ref. [25] that as u — > the normal phase conductivity satisfies^] 

Ac? r_i_ 12r 4 — o 2 

lm(a n )= / -± + 0(1) where T= + { . (4.20) 
3(4r_J. + p A ) uj 167r? + 

From our numerics, we can study this pole for T < T c . We find that there is an additional 

contribution to the pole below the critical temperature. Close to T c , it takes the form 

Im(a) = ImK) + -(T c - T) . (4.21) 

The coefficient A depends on the charge of the scalar field and its boundary condition at 
infinity. Consider, for example, the case where q = 3. For the dimension one operator, 
A ~ 15 and for the dimension two operator, A ~ 12. Therefore d!m(a)/dT is discontinuous 
across T c . In computing these slopes, p is held fixed as T is varied. 

Given that we appear to see an exact gap emerge as T — > 0, i.e. that Re(cr) vanishes 
identically for u> < lo 9 at T = 0, we can expect that at small but finite temperatures, the 



9 This pole vanishes in the probe limit since that limit requires pocl/gasg^oo. 
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zero frequency limit of the real part of the conductivity is governed by thermal fluctuations 

lim Re(cr) ~ e~ At/T (4.22) 

when Aj/T 3> 1. We can further ask whether Aj is related to uj g . If we parametrise this 
relation by A, = aiU g then, numerically using the lowest temperatures accessible to us, we 
find that 



3 
6 
12 



0.13 0.14 

(4.23) 

0.34 0.30 



0.45 0.44 

From ref. [3], we know that in the probe limit, which corresponds to large q, these values of 
at should become close to 1/2. Extracting these values numerically is delicate as one needs 
to obtain very low temperatures. The accuracy of our numerics decreases with q and we 
have not quoted the q = 1 values because our numerical results are insufficiently robust. 



On the other hand, the exponential behaviour (4.22) is clearly seen, as is the fact that the 
ai are less than 1/2 for finite q. 

4.3 Comments on gaps and pairing 

The fact that the conductivity in the superconducting phase tends toward the value of the 
normal phase conductivity as to —* oo indicates that the degrees of freedom responsible for 
high frequency conductivity are those of the normal phase. 

At zero temperature and frequencies u < uj g conduction is non-dissipative, i.e. the real 
part of the conductivity vanishes identically, at least within the accuracy of our numer- 
The region where the conductivity vanishes implies a corresponding gap in the charge 
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spectrum — if there were asymptotic charged states with energy below u g , they would 
contribute to the conductivity at these low frequencies. Yet this statement is not entirely 
true: we know that there is a Goldstone boson in our symmetry broken phase. Usually 
a Goldstone boson leads to a nonzero conductivity all the way down to cu = due to the 
presence of multi-Goldstone boson states. The fact that this does not occur in our models 
is possibly a large N effect. 

In a standard weak coupling picture of superconductivity, the gap uj g is understood as 
the energy required to break a Cooper pair into its constitutive electrons. The energy of 
the constituent quasiparticles is given by A and then uj g would be some integer multiple of 



"'Recall we defined uj g to be q{Oi) or yj q(02), motivated by figure Q 
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this energy. The non-integer relation between A and u g in table (4.23) shows that we are 
clearly not in a weak coupling regime and that such a quasiparticle picture is not applicable, 
except perhaps in the q — * oo limit, in which we recover the probe result |3] oj g = 2A. At 
strong coupling we would expect to be closer to a Bose-Einstein condensate scenario than 
to a BCS like weak pairing description. On the other hand, a non-integer relation is almost 
inevitable if tu 9 /A is to depend continuously on the charge q. It may be that microscopic 
realisations of holographic superconductors, that is, embed dings of our setup into string 
theory, will place constraints on the masses and charges of the scalars that condense. 

It is interesting to note that for the dimension two condensate and charges q > 3, even 
though a.i changes by at least a factor of three, uJ g /T c remains close to the value 8. The 
constancy of co g /T c can be seen by comparing Fig. 4b and Fig. lb. We conclude that the 
energy to 'break apart' the condensate is insensitive to the charge for q > 3. In other words, 
the probe limit is reached rapidly for this ratio. 



5 Critical magnetic fields 

In this section, we start our investigation of the effect of magnetic fields on our holographic 
superconductor. We will argue that our model behaves as a type II superconductor. Recall 
that the difference between type I and type II superconductors lies in the way the Meissner 
effect disappears as the temperature of the material is raised, as we now review. 

The Meissner effect is the observation that at low temperatures superconductors expel 
magnetic field lines. The existence of this effect implies that there is a critical magnetic 
field H c above which the superconducting order is destroyed and the material reverts to its 
normal state: The superconductor must perform an amount of work H 2 V^/8ir to expel an 
applied field H from a volume V3. The critical field strength is then obtained by equating 
this work to the difference in free energies between the normal and superconducting states 
of the material. For H > H c it is no longer thermodynamically favorable to superconduct. 

How the superconductor reverts to a normal state as the magnetic field is increased 
depends on the nature of the material. For type I superconductors, there is a first order 
phase transition at H = H c , above which magnetic field lines penetrate uniformly and the 
material no longer superconducts. For type II superconductors, vortices start to form at 
H = H c \. In the vortex core, the material reverts to its normal state and magnetic field lines 
are allowed to penetrate. The vortices become more dense as the magnetic field is increased, 
and at an upper critical field strength H = H C 2, the material ceases to superconduct. 
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In determining whether our model is type I or type II, we are faced with the limitation 
that currents in the model do not source electromagnetic fields. One immediate consequence 
of this limitation is notational. As the material does not produce its own magnetic fields, 
the applied magnetic field is the actual magnetic field and we can set H = B. Given the 
limitation, the Meissner effect cannot strictly speaking exist: to exclude the magnetic field, 
the current produced by the external magnetic field must produce an equal and opposite 
canceling field inside the sample. However, we will show in a later section that holographic 
superconductors do generate the currents required to expel magnetic fields (the London 
equation) and that the theory can consistently be weakly gauged. Therefore we assume for 
the moment that (a gauged version of) our model attempts to expel fields in the usual way 
for superconductors. 

This brings us to a subtlety. We would like to work with a 2+1 dimensional model 
interacting with a 3+1 dimensional electromagnetic field. This is a realistic setup for a thin 
film superconductor. We apply the magnetic field normal to the material. Assume for the 
moment that the 2+1 dimensional sample is a disk of radius R. In order for the disk to 
expel the magnetic field, the disk must produce a current circulating around the perimeter. 
Solving Maxwell's equations, this current will expel a field not only in the area nR 2 of the 
disk but in a larger volume of size V3 ~ i? 3 . As mentioned before, the amount of work that 
the superconductor must do to exclude an applied magnetic field from a volume V3 scales as 
H 2 Vs. In the large R (thermodynamic) limit, the superconductor does not have enough free 
energy available to expel a magnetic field from such a large region; the difference in the free 
energies between the normal and superconducting phases can scale extensively only as R 2 . 
Thus magnetic fields of any non-vanishing strength will penetrate a thin superconducting 
film and B c = 0. This argument is illustrated in figure [5] below. 

So far, our discussion applies to any thin superconductor in a perpendicular magnetic 
field. For our model to be type II, we still need to establish that the material remains a 
superconductor for values of the magnetic field up to some upper critical field strength B C 2. 
To that end, we now investigate the largest value of the magnetic field for which our scalar 
condenses. 

5.1 Superconducting droplets 

In this section, we will examine the effect of a constant background magnetic field on the 
dynamical ability of the condensate to formj^j] For each temperature T < T c , we expect 
llr The computations in this subsection are similar to those in [7], but our interpretation is different. 
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Figure 5: In order to prevent the flux from penetrating the superconductor, of area i? 2 , 
the currents would have to do enough work to expel the field from a volume of size i? 3 , 
as shown in the left hand figure. This work cannot be supplied by the free energy gain of 
superconducting on the thin film. Therefore, the flux always penetrates the film, as shown 
in the right hand figure. 

there to be some critical field strength B C 2 above which the condensate cannot form. We 
find B C 2 by starting in a phase with a large magnetic field and no condensate, and showing 
that this phase develops an instability towards condensation of the scalar as B is lowered. 

Our starting point is a dyonic black hole background. This is a Reissner-Nordstrom 
AdS black hole with both electric and magnetic charges and no scalar hair. As we have 
already seen, the electric charge of the black hole gives the charge density of the field theory. 
The magnetic charge gives the value of the background magnetic field, as explained in, for 
instance [20J. The form of the solution is well known (see e.g. [20]). The metric takes the 



form (3.1 ) with x = and 



g(r) = r 2 - — (44 + P 2 + B 2 ) + ^(p 2 + B 2 ) . (5.1) 
In this section, we find it convenient to work in polar coordinates dx 2 + dy 2 = du 2 + u 2 dip 2 . 



To the vector potential (3.2) we have to add a magnetic component: 

A = p ( — - - ) dt+ -Bu 2 d^. (5.2) 
\r+ r) 2 

The horizon radius r + is determined implicitly by the temperature via 

12ri -p 2 -B 2 

T= ± — —„ . 5.3 

16vrri v ; 
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We will choose a gauge in which the scalar tp is real. Given the above background, we will 
treat the scalar field as a perturbation and look for a solution that is well behaved at the 
horizon (no logarithmic divergence) and chooses one of the two fall-offs, = or ^ = 
at the boundary, corresponding to our two choices of boundary operator. We interpret the 
existence of such a zero mode solution as the onset of an instability for the condensate of 
the corresponding scalar operator to form. 

We are interested in static, axisymmetric solutions in which all fields are independent 
of t and 93. The scalar field thus satisfies the equation 



1 



d u (ud u ip) + d r (r 2 gd r ip) + 



q 2 p 2 

2 



(r — r + Y 



,2, ,2 



-B + 2r z 



= 



(5.4) 



This linear PDE can be solved by separation of variables: 



ip(r,u) = R(r)U(u) 



(5.5) 



where U{u) solves the equation for a two dimensional harmonic oscillator with frequency 
determined by B: 

U'' + \u'-(^fju = -XU, (5.6) 

R satisfies 



(r 2 gR'Y + 



■" 2 p 2 (r — r + ) 2 r, 



R — XR , 



(5.7) 



and the separation constant A = qnB. Clearly, the condensate is now clumped, with a 
Gaussian profile. We expect that the lowest mode n = 1 will be the first to condense and 
lead to the most stable solution after condensing. Therefore we choose 



U (u) = exp(— qBu 2 /4) 



(5.8) 



From the equation of motion (5.7), it is straighforward to analyze the near boundary 



and near horizon asymptotics of R{r). The near boundary behavior is the same as was 



found above in (3.11). The near horizon behavior on the other hand takes the form 



R = cq + ci ln(r/r^ 



1) + ... 



(5.9) 



We are looking for a regular zero mode and thus set c\ = 0. The equation for R is linear 
and so cq is arbitrary and can be set to one. Requiring either ip^ = or t// 2 ) = at the 
boundary thus produces a curve of solutions in the (p,B) plane. We call this curve B C 2{p) 
and have plotted it in figure |6j There will be a localised droplet of condensate on the lower 
left region of these figures. 
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Figure 6: We plot the critical field B C 2, below which a droplet of condensate forms, versus 
temperature: a) The 0\ case with, from right to left, q = 12, 6, 3, and 1; b) The 02 case 
with, from right to left, q = 12, 6, 3 and 1. In the lower left region there are circular droplets 
of superconducting condensate. In the top right region there is no superconductivity. 

It is important to note that in order to see the onset of this instability we only needed 
to work to first order in the condensate, which is small just below the transition. We will 
see in a later section that at next order the condensate causes magnetisation currents. A 
proper treatment of these currents would have to include the effect of the backreaction 
of these currents on the external magnetic field via Maxwell's equations. For instance, at 
sufficiently low temperatures the superconducting droplets presumably grow and trap the 
magnetic flux into vortices. Magnetic screening is an important feature of vortex physics, 
yet this would take us beyond the AdS/CFT model. It is fortunate, therefore, that these 
effects are not important at the transition itself. 

5.2 Thermodynamics 

In this section we investigate quantitatively the effect of the magnetic field on the free 
energy of our model. We are not able to carry out the most obvious calculation, which is a 
determination of the magnetic field dependence of the free energy in the superconducting 
phase, since we don't have the general solution for a hairy black hole with both electric and 
magnetic charge. However, we will see clear evidence that the transition is second order at 
B = 0, where we can compute the free energy. In passing, we also remark on the strong 
diamagnetism of our material in the normal phase. 

We choose to work in the canonical ensemble, at fixed charge density. Recall the ther- 
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modynamic identity 

E + PV = ST + pQ, (5.10) 

where E is the total energy, P is the pressure, V the volume, S the entropy, and Q the total 
charge. The combination — PV is also the value of the potential function Q in the grand 
canonical ensemble. We would like to work instead with the free energy F in the canonical 
ensemble 

F = n + pQ = -PV + pQ = E - ST . (5.11) 

For our particular theory, we have an additional relation that comes from the tracelessness 
of the stress-tensor. However, we have to be a little careful here because in the presence 
of a magnetic field B, it is possible to define two different pressures. The diagonal spatial 
components of the stress-tensor are related to P via a magnetization M = mV, T n = 
P — mB. Tracelessness thus implies E = 2{PV — MB), and we may write 

F(B,0) = ~ + y.Q-MB . (5.12) 

First, we would like to compute the value of the free energy of a configuration with a 
magnetic field but no condensate. We wrote the corresponding gravity background, a dyonic 



black hole, in (5.1) above and the temperature in (5.3). To determine fi, one calculates 
the on- shell value of the regulated Euclidean action. The action is regulated by the usual 
Gibbons-Hawking term and a boundary cosmological constant. We do not reproduce the 
details here as they can be found in [20]. The result is 

n(B,0) = ^(-r\- P ^ + ^]v . (5.13) 



r+ V 4 4 

To determine F, we must add pQ = p 2 V: 



F(B,0) = ^-(-rl + 3 l + 3 ^)V. (5.14) 



The magnetic field dependence of (5.14) is unusual and means that, while not super- 
conducting, the normal phase of the material is strongly diamagnetic at low temperatures. 
Note that at small values of p and B, we may replace r + with 47rT/3. Thus the mag- 
netic susceptibility \ = d 2 F/dB 2 becomes of order 1/T. The dimensionless quantity yT 
is naively of order one for this model. If we normalize the action to be consistent with 
established AdS/CFT dualities, such as the M2-brane theory, then the susceptibility will 
scale with a power of N in the large N limit. Compare this result with a typical 3+1 di- 
mensional, non-ferromagnetic metal. In 3+1 dimensions, the susceptibility is dimensionless, 
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and approximating a metal as a free electron gas, the susceptibility is suppressed by a power 
of the fine structure constant and is typically tiny. There is a general lesson here: quantum 
critical theories will often be strongly magnetic because there is no small coupling or scale 
to suppress the magnetic susceptibility. 

The smallness of the susceptibility for a 3+1 dimensional electron gas explains an ap- 
proximation that is typically made in calculating H c for superconductors. Namely, the 
dependence of the free energy of the normal phase on the magnetic field is neglected. Be- 
cause of the large diamagnetism of our model, we clearly would not be able to make this 
approximation. Nevertheless, the argument above that B c = for our model still holds. 
The reason, as we explained above, is that to be a perfect diamagnet, the superconductor in 
2+1 dimensions essentially has to have an infinite susceptibility that scales with the system 
size R. 

Next we compute the free energy of the system with a condensate and no magnetic field. 



The parameter e in the expansion (3.19) and (3.20) is naturally interpreted as the energy 
density in the boundary field theory, eV = E. Meanwhile, we can obtain p, and Q = pV 



from the asymptotic expansion of (f> (3.10). Putting the pieces together, we have from (5.12 ) 

F(0,O)=(~ + w)v . (5.15) 



2 

This result agrees with the value of the regulated on-shell Euclidean action calculated in 



(3.31 ). The value of the regulated action is f2 = — eV/2. 

Being careful to compare (5.14) and (5.15) at fixed T and p, figure [7] displays the free 
energies for representative values of the parameters. Note the continuous second order phase 
transition between the normal and superconducting phases at B = 0. In order to show the 
continuity of the transition at finite B, we would need the black hole background with both 
condensate and magnetic field, at least to second order near the critical temperature. The 
qualitative similarity between the instability with and without magnetic fields suggests that 
the transition will be second order in general. 

Having computed the free energy, it is straightforward to obtain the specific heat by 
differentiating: c = —T/Vd 2 F/dT 2 . At low temperatures we do not find the exponential 
suppression of the specific heat typical of s wave superconductors. Rather the specific heat 
vanishes as a power law as T — > 0. It is difficult to determine the precise power due to 
numerical sensitivity at low temperatures. A likely source of this power law behavior is the 
presence of a Goldstone mode in our system. 

A power law rather than exponential behavior at low temperatures is also observed 
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Figure 7: The solid blue line is the free energy of the hairy black hole {B = 0). The 
dotted red line tangent to this blue line is the free energy of the electric black hole with no 
condensate {B = 0). The dotted red line at the top is the free energy of a dyonic black hole 
(no condensate) with B/T 2 = 20. The plot is for 0\ with q = 3. We see that the magnetic 
field raises the free energy of the normal phase. 

for the quantity n s {T) — n s (0). Here we define n s to be the coefficient of the pole in the 
imaginary part of the conductivity as uj — > 0. 

6 Magnetically induced currents in the superconducting phase 

In this section we continue our study of the superconducting phase in the presence of a 
finite magnetic field. The key physics we wish to examine are the currents generated by the 
background magnetic field. These currents are responsible for the Meissner effect once the 
theory is coupled to dynamical photons. 

There will be two main discussions in this section. Firstly, we shall exhibit the London 
equation analytically at low temperatures. Secondly, we study the phase diagram of the 
(ungauged) theory in the presence of a homogeneous background magnetic field. 
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6.1 The London equation at low temperatures 



We would like to explain how the London equation arises in our model The London 



equatior 13 



Ji(u, k) = —n s Ai(u>, k) (6.1) 

was proposed (in a gauge where the order parameter is real) to explain both the infinite 
conductivity and the Meissner effect of superconductors. This equation is understood to 
be valid where u and k are small compared to the scale at which the system loses its 
superconductivity. In our case, that scale will be (Oi). One important and subtle issue in 
understanding this equation is that the two limits u — > and k — > do not always commute. 
In the limit k = and u) — > 0, we can take a time derivative of both sides to find 

Ji(w,0) = ~Ei(u,0) (6.2) 

LO 

explaining the infinite DC conductivity observed in superconductors. On the other hand, 
in the limit u = and k — > 0, we can instead consider the curl of the London equation, 
yielding 

ie i:j ik j J l (0,k) = -n,Bi(p,k) . (6.3) 

Together with Maxwell's equation e^ l djBi = Air J 1 , this other limit of the London equation 
implies that magnetic field lines are excluded from superconductors. 

Thus far in the paper, we have explored the first limit, having set k = and explored 
the frequency dependence of the conductivity. We would now like to argue that the London 
equation holds more generally, including in the limit where u is sent to zero first. To make 
life easier, in this section we shall work in the probe limit (q — > oo) in which the metric is 
kept fixed to be simply the Schwarzschild AdS black hole. In the probe limit, the scalar and 
Maxwell field form a decoupled Abelian-Higgs system in this background. As mentioned 
previously in section |4j by decoupling the metric fluctuations, we will remove the additional 
divergence in the conductivity at u — > due to translation invariance. The background 
metric is 

dv 2 

ds 2 = -g(r)dt 2 + — + r 2 (dx 2 + dy 2 ) , (6.4) 

g(r) 

where g(r) = r 2 — r\jr. 

Assume that we have solved self-consistently for At and ip in this Schwarzschild back- 
ground. We then allow for perturbations in A x that have both momentum and frequency 

12 There is some overlap of this section with [S] which appeared as we were completing this work. 
13 Recall that we have defined n s to be the coefficient of the pole at uj = in lm(cr). 
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dependence of the form A x ~ e - lujt + tk y. We have taken the momentum in a direction or- 
thogonal to A^. This allows us to consistently perturb the gauge field without sourcing any 
other fields. With these assumptions, the differential equation for A x reduces to 

- - ^) A x + (gA'J = 2q^ 2 A x , (6.5) 
9 H J 

where ' denotes differentiation with respect to r. Ignoring the radial dependence, this 
equation describes a vector field with mass proportional to q 2 ip 2 . This mass, which is 
symptomatic of an underlying bulk Higgs mechanism, should give rise to the usual effects 
of superconductivity, but the radial dependence and the AdS/CFT dictionary cloud the 
intuition. 



To clarify the situation, first note that we have already explored this equation (6.5) 
numerically in the case k = both here in section [4] and also previously in [3] . In section 
[4j strictly speaking we were not working in the probe limit, but we found that as q — > oo, 
the results approached those of [3] where we were indeed working in the probe limit. The 



results indicated the imaginary part of the conductivity had the form (6.2), verifying the 
London equation in the limit k = and u — > 0. 



The next observation is that given the structure of (6.5), the limits u> — > and k — > 



must commute. To compute n s , we simply set both to and k to zero and solve (6.5). Thus 



in the probe limit we directly obtain the magnetic London equation (6.3). Away from the 
probe limit, these limits may fail to commute. As we saw in section El Maxwell's equation 



for A x (4.1) depends also on the metric fluctuations gt x . We saw that for k = and 
lo 7^ 0, the metric fluctuations could be replaced by an additional effective mass term for 
A x . However, for to = and k ^ 0, more metric and gauge field fluctuations were sourced, 
and the differential equation governing gt x are no longer first order. 

We now attempt to approximate n s analytically at very low temperature, following a 
method that was used successfully in [3J. Without knowing an analytic form for ■0, it is not 
possible to provide an exact solution to this differential equation for A x . Nevertheless, we 
find numerically that in the dimension one case ^/2tp ~ {0\)/r to a good approximation 
everywhere. At very low temperature, r + — > and our background approaches AdS in 



Poincare coordinates. Introducing a new radial variable z = 1/r, eq. (6.5) reduces to the 
Klein-Gordon equation with mass proportional to (Oi): 

{uj 2 -k 2 -q 2 {0 1 ) 2 )A x + A x = . (6.6) 

Here a dot denotes differentiation with respect to z. We are implicitly working at low 
frequencies where u> 2 , k 2 <C q 2 (0\) 2 . Since the horizon is at large z, we impose the boundary 



32 



condition that A x be well behaved there to find 

A x = a x e- iut+iky ~ Xz , (6.7) 

where A 2 = g 2 (0i> 2 + k 2 - uj 2 « g 2 (Oi) 2 . 

To obtain the conductivity we expand A x near the boundary z = in the low frequency 
case: 

A x = a x {l - \z + 0{z 2 )) . (6.8) 

From the AdS/CFT dictionary, described several times above, we can interpret the zeroth 
order term as an external field strength and the linear term as a current J x . Thus, this 
expansion gives us a modified London equation: 

J x = -vV(0l> 2 + fc 2 -^ 2 a x . (6.9) 

In the limit u>, k <C q(Oi), we get precisely the London equation: 

J x = -q{0 1 )a x . (6.10) 

We have verified numerically that the strength of the pole in the imaginary part of the 
conductivity is indeed very close to q(0\) at low temperatures. 

A similar estimate of n s for the {O2) theory is given in Appendix B. Since the approxi- 
mation \/2i/) ~ (O2) /r 2 is not as good, the estimate is off by about 25%. But the main point 
is that by allowing for both a momentum and frequency dependence, we are free to choose 
the order of limits in which we send the frequency and momentum to zero. If we set u> = 
first, then we are manifestly dealing with a purely magnetic external field. The London 



equation (6.10) is precisely what we need to describe both the expulsion of a magnetic field 
from the superconductor and the infinite conductivity. 

The London equation leads to the magnetic penetration depth 

A2 = i- < 6 - n > 

via the Maxwell equation for the curl of the magnetic field: 

- \7 2 B = V x (V x B) = 4vrV x J = -4vrn s V x A = -4irn s B . (6.12) 

Therefore 

V 2 B = -^B, (6.13) 

implying that static magnetic fields can penetrate a distance A into the superconductor. 
Although we argued above that a 2+1 dimensional superconductor cannot expel a per- 
pendicular magnetic field, this lengthscale will still play an important role in a gauged 
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extension of our model. For instance, at low temperatures we expect the flux to be confined 
to vortices, and their size will be determined by A. Some comments on weakly gauging a 
holographic superconductor appear in section 7 below. 



6.2 The superconductor in a finite magnetic field 

In the previous section, we studied the onset of superconductivity in the presence of a 
constant background magnetic field. We found that at sufficiently low temperature there 
is a second order transition to a superconducting droplet. Thus the first consequence of 
the magnetic field is to confine the superconducting condensate to a finite region. In this 
section we are interested in characterising the currents associated with the droplet phase. 



The formalism here is similar to that used in section 5.1 Let us write the background 
metric in polar boundary coordinates: 

ds 2 = -g(r)dt 2 + d ^ + r 2 (du 2 + u 2 d^ 2 ) . (6.14) 

g{r) 

We will choose a gauge in which the scalar ip is real. We are interested in static, axisym- 
metric solutions in which all fields are independent of t and (p. In this case, it is consistent 
to set A r = A u = 0. The Maxwell and scalar field equations become 

-d u (ud u A t ) + gd r (r 2 d r A t ) = 2r 2 q 2 ^ 2 A t , (6.15) 
u 

ud u (-d u A v ) + ^drigdrA^) = 2r 2 q 2 iP 2 A v , (6.16) 
u 



-d u (ud u il)) + d r (r 2 gd r ^) + 
u 



ip = 0. (6.17) 

g ' w r J 

To study the complete effects of a magnetic field on the superconductor, one would have to 
solve these coupled nonlinear partial differential equations. This will require more sophisti- 
cated numerical methods. 

We will now study these equations in two different limits in which they become tractable. 
Firstly we look at the case where the magnetic field is small compared to the background 
charge density. In this region we can treat the magnetic field as a perturbation and linearise. 
Secondly, we look at temperatures just below the formation of the superconducting droplet. 
Here we can treat the scalar field as a perturbation. Finally we will put together the main 
features of these two limits to obtain a qualitative picture of the full phase diagram. 
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6.2.1 Small magnetic fields 



Consider the case that the magnetic field is weak and can be treated as a perturbation of 
the solution with no magnetic field. We shall see that a static magnetic field generates 
currents in the superconductor. 

Assume that we have solved self-consistently for At and ip in the absence of a magnetic 
field, as we have been doing in earlier sections of this paper. We introduce a small but 
nonzero in this background. The equation (6.16) can now be separated A v = uV{u)S{r) 
with 



" 1 , 
V H V + | (- 

u 



1 



V = , 



(gs' 



% + 2gV ) S = , 



(6.18) 



(6.19) 



where c is a constant of integration. The solution for V which is regular at the origin is just 
the Bessel function J\{cu) where c is real. Note that c can be chosen at will and its value 
determines how spread out the magnetic field is on the boundary. As c — > 0, there is an 
arbitrarily large region centered at the origin where A v oc u 2 , corresponding to a uniform 
magnetic field. 

We now need to solve for S in order to find the current. One cannot solve for S 
analytically, but as usual, at large r, S = + S^/r where SW are constants and 



<?(1) 

J v = uJ 1 {cu)S {1) = —A, 



(6.20) 



S(o) v ' 

In the absence of a condensate, it is easy to see that the current vanishes. Subtracting a 



constant from (6.19) we can view this as an equation for S = S — which vanishes at 



infinity. If ip = 0, we can multiply (6.19) by S and integrate over the region outside the 



horizon. The result is a nonpositive integrand which integrates to zero, implying S = 0. 



The equation (6.19) is identical to the Maxwell equation (6.5) we wrote previously, with 
lo = and k 2 — ► c 2 . In fact, the only difference between this calculation and the one in 
section (6.1) is that previously we used cartesian coordinates and a standard plane wave 
dependence for the magnetic field, while here we used polar coordinates and the associated 
Bessel function dependence. Therefore at low temperature, the coefficient of the London 



equation (6.20) in the limit c — * will be identical to that of the previous subsection. Of 



course, since we are not trying to solve for S analytically, we can now consider all T < T c , 
whereas previously we were restricting attention low temperatures. 
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6.2.2 Small condensates 



We now revisit the calculation in section 5.1 where we examined the limit where ip is very 



small, i.e. we are near T = T c . The right hand sides of (6.15) and (6.16) vanish, so 



satisfies the source-free Maxwell equation. One solution to (6.16) is then simply 

_Bu 2 



(6.21) 

corresponding to a uniform magnetic field. Similarly, At can be just a function of r. As 
before, the equation for tp(r,u) = R{r)U{u) can be separated, yielding the differential 



equations (5.6) and (5.7). The solution for U is a Gaussian profile as in (5.8). Even though 



the condensate is nonzero, there is no current generated in this limit because we are working 
to first order in i/j and A^ is independent of r. 

We now wish to support our general picture of currents induced by the external magnetic 
field by extending our small tp expansion to C(V ;2 ) m order to compute the induced current. 
Note that unlike in section 5 above, we are in the probe limit here. Substituting ( 6.21|5.5 ) 



into the right hand side of (6.16) we get 



ud u {-d u 8A v ) +r 2 d r {gd r 5A tp ) = q 2 Bu 2 e- qBu2/2 r 2 R 2 (r 



(6.22) 



This equation doesn't separate, but we believe that we can understand its solutions as 
follows. Immediately we can see that there is a source term for A v that will cause an r 
dependence in the solution and hence produce a current in the field theory. Note that 
regularity near u = requires that for small u 



5A V = u 2 S(r) 



(6.23) 



This already shows that the current (6.20) grows with radius u away from the center of the 



condensate. To understand the structure of the solution at large u, we try 

5A V = u p e- qBu2 / 2 S p (r) . 



(6.24) 



The first term in (6.22) produces terms with u dependence given by the exponential times 
u p+2 , u p and u p ~ 2 . Clearly p cannot be positive since there are no terms to cancel the 
highest power of u. Setting p = leads to the approximate large u solution 

1 



SA = _ e ~iBuV2 r 2 R 2 {r) 
B 



This solution can be systematically corrected by taking 

oo 

SAp = £ n- 2m e-^ 2 / 2 S 2m (r) . 



(6.25) 



(6.26) 



m=0 
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The functions S2 m are determined iteratively and algebraically in terms of the lower m 
functions and their derivatives. This shows that the current dies off exponentially far from 
the condensate, as one expects, since persistent currents must vanish when there is no 
condensate. To make this discussion more rigorous, one should show that there is indeed a 
solution to the equations matching these two asymptotic behaviours. 

6.2.3 Comments on the full phase diagram 

By piecing together the results from the two limits we have just discussed, small magnetic 
field and small condensate, we arrive at the following picture for the superconductor in a 
constant external magnetic field, i.e. without dynamical photons. Because we have not 
solved the full equations, the description that follows is a minimal interpolation between 
the different regions we have studied. This phase diagram is sketched in figure [8] below, 
with the regions that are accessible to a linearised analysis shaded. 




B 

Figure 8: Schematic illustration of the full phase diagram of the theory at finite temperature 
and external magnetic field. The shaded regions indicate where we can obtain a description 
of the phase by linearising the full equations in either the magnetic field or the condensate. 
In the low temperature, low magnetic field phase, the condensate is concentrated in a 
circular droplet with circular currents in the superconducting region. 

Firstly, for any finite magnetic field the superconducting condensate will be localised 
to a finite circular region. As the magnetic field becomes smaller the region grows until 
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it occupies the whole plane in the B — > limit. In general there are exponential tails of 
superconductivity reaching out to infinity. It would be interesting to ascertain whether in 
the zero temperature limit the condensate becomes completely localised. 

Secondly, for any finite superfluid density and magnetic field there are always circular 
currents generated. The currents are largest just inside the boundary of the condensate and 
die off exponentially at large radii. When coupled to dynamical photons, these currents act 
to expel the applied magnetic field. 

7 Photon mass, symmetry breaking and infinite conductivity 

In this section we will make some more formal observations about holographic supercon- 
ductors. In particular, we wish to give an interpretation of the pole in the conductivity at 
ijj = in the superconducting phase. By explicitly coupling the theory to a photon, we will 
show that the pole is directly related to the photon becoming massive, 14 

As we have mentioned several times by now, the 2+1 dimensional theory we have been 
considering does not have a dynamical photon. The £7(1) symmetry that is spontaneously 
broken is global. Although the dynamics of spontaneous symmetry breaking does not de- 
pend on the photon, much of the interesting phenomenology of superconductors is concerned 
with the interaction of the theory with a dynamical photon. 

The 2+1 theory can be coupled to a photon through the standard J^A 11 interaction. 
To make the photon dynamical, we can add an F 2 term to the action, with F = dA. 
Electromagnetic phenomena such as screening are determined by the effective action for 
the photon. We can obtain this action by integrating out all the other degrees of freedom. 
In terms of the Euclidean partition function, we have 



In these expressions X denotes the degrees of freedom in the 2+1 dimensional theory. 

14 There is a temptation to identify the pole at uj — in the conductivity as being due to the Goldstone 
boson of the spontaneously broken (global) (7(1) symmetry. This is not correct however. The Goldstone 
boson is manifested as a pole in the retarded Greens function of the current at uj 2 — k 2 — 0. Whether this 
pole persists at uj = k — depends on the direction in which zero is approached in the (ui, k) plane. We are 
setting k = first. In this order of limits there is not a pole at uj = 0. This is manifest in our results: if there 
had been a pole in the retarded Greens function at uj = 0, then the conductivity would have had a double pole 
at u — 0, because of the relation <j(uj) = —iG^j w j m (uj)/uj [From Ohm's law: J x — aE x — iujaA x — Gj j x A x -\ 



Z 




(7.1) 



(7.2) 
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Unlike in previous sections, in this discussion we are taking the photons to also be 2+1 
dimensional. This is because we are not interested in mimicking experimental setups, but 
rather in demonstrating a formal property of the theory. Note that the coupling e 2 is 
therefore dimensionful. 

Up to quadratic order in the Maxwell field, the effective action can be straightfor- 



wardly obtained by expanding out the exponent in (7.1), integrating over X, and then 
re-exponentiating. For our theory, in which there is a charge density but no background 
currents, we get 

Ses.[A] = ^2 f d 3 xF^F^ + l -f d 3 xd 3 y(J»r)c(x - y)A fl (x)A„(y) + JdtQfi. (7.3) 

This expectation value is evaluated in the theory without dynamical photons. The last term 
arises from setting 

d 3 xpA t (x) = dtQfi. (7.4) 



This term is an overall constant, does not affect the dynamics, and so we shall drop it from 
this point on. 

In expanding the effective action in powers of A one might worry about the fact that we 
have dropped the higher order, interacting terms. We are helped here by the large N limit, 
which is a classical limit of the theory. The coefficients of the higher order terms in A in the 
action are given by connected higher n-point functions of the current: {J n ) c - Connected 



diagrams arise upon re-exponentiating the partition function (7.1 ) after integrating out the 
CFT degrees of freedom. These connected diagrams are computed in AdS/CFT by using 
Witten diagrams in the bulk theory which have vertices or loops or both. Therefore they 
are suppressed by inverse powers of N compared to the disconnected diagrams at the same 
order. If we were to rescale A so that the coefficient of the A 2 term in the action were 
Oil), then the higher order terms in the effective action would be suppressed by powers 
of N. In this paper we are studying the theory as a function of frequency at zero spatial 
momentum. For the F 2 term in the action also to be bigger than these interaction terms 
we need to 2 /e 2 > j d 3 x(JJ) c . 



An immediate feature of the effective action (7.3) is the possible presence of a mass 
term for the photons, via (J^J U ) C . Let us see if we can extract the photon mass. Lorentz 
invariance is broken by the charge density, so we need to define what we mean by the photon 
mass. From our results we can obtain the energy of photons in a frame where they are at 
rest relative to the background charge density. It is reasonable to associate this energy with 
the mass of the photons. Therefore, to find the photon mass we need to exhibit an on shell 
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photon mode with k = 0. The energy uj of this mode will be the photon mass, m 7 . An 
example of a physical consequence of this definition of the photon mass is that generic zero 
momentum processes involving massive photons will decay in time like e~ m ~' t . 

Because we have a quadratic effective action, the spectrum can be obtained directly from 
the classical equations of motion. It is straightforward to obtain the equations of motion 



that follow from the effective action (7.3). Let us note that it is consistent to restrict to a 



mode in which A x is the only nonvanishing component and where 

A x (uj) = yB^ 1 ■ (7-5) 

The y dependence has been included so that the mode is not pure gauge at u = 0. Rather, 
it reduces to a constant magnetic field. We then Wick rotate so that we are in Lorentzian 
signature. The equation of motion for this mode is 

{uj 2 + e 2 G% JX {oj)) A x {uo) = 0. (7.6) 

Here we used the fact that the retarded Green's function in momentum space is the Fourier 
transform of the Euclidean Green's function in momentum space, up to possible contact 
terms. However, it was shown in [25] that the contact terms are not present in the case of 
the current-current correlator in this theory. We are also assuming that G^jxjy — 0. 
Using Gj- xjx (uj) = iua{uJ), we can recast the 'dispersion relation' as 

uj (uj + ie 2 o(uj)) =0. (7.7) 

The solution to this equation gives the photon mass m 7 = u. We can see immediately that 
uj = is a solution provided that the conductivity does not have a pole (or worse) as uj — > 0. 
In particular, if the conductivity tends to a constant or vanishes in the uj — > limit, then 
the photon is massless. In Lorentz invariant theories we expect the photon to be massless 
in phases where the electromagnetic U(l) is unbroken. Conversely, we expect the photon to 
gain a mass if the electromagnetic symmetry is broken. If the imaginary part of a(uS) has a 
pole as uj — > 0, the photon will have a nonzero mass. It is important to note, however, that 
we also have a medium, the charge density, which breaks Lorentz invariance. 

How do these expectations compare with our results? Take the superconducting phase 
first. There indeed the imaginary part of the conductivity has a pole at the origin, and so 



uj = is not a solution to the dispersion relation (7.7), and hence the photon is massive, 
perfectly consistent with expectations. We found above that o~{uj) ~ in s /uj as uj — > (this is 
the definition of n s for us). This behaviour is not exact. At zero temperature, for instance, 
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the imaginary part of the conductivity goes to zero at the gap u = u g . In general we will 



have to solve (7.7) numerically. However, if e is small enough we will have 

m-y = e^JrT s provided e^/rig <C 0J g ■ (7i 



In the probe limit (recall this was q — * oo), the conductivity a is constant in the normal 
phase. There is no pole as uj — > and therefore we find that the photon is massless, again 
consistent with our naive expectations for a symmetric phase. 

Beyond the probe limit, when back reaction of the scalar field on the metric is taken 
into account we found that the conductivity also diverged as ui —* above T c . This pole 
in the conductivity was not due to spontaneous symmetry breaking, but rather due to a 
translationally invariant charged medium and can be moved away from uj = by impuri- 
ties or by considering finite momentum. A charged medium often results in screening of 
electromagnetic fields, and our pole at finite u is consistent with this. We should also note 
that the limits uj — ► and B — > do not commute when back reaction is included |25| |6] 
and so it is subtle to take simultaneously the small uj and B limits in this case. 



8 Comparison to Landau-Ginzburg 

Before concluding we would like to discuss the microscopic status of our phenomenological 
holographic superconductor, in particular the extent to which it is similar and distinct from 
a dressed-up version of Landau-Ginzburg theory. At a first glance what we are doing feels a 
lot like Landau-Ginzburg theory — we have studied a theory like the Abelian-Higgs model, 
albeit in one higher spacetime dimension and with a dynamical black hole metric, in which 
there is a condensate for a charged complex scalar field. Furthermore, in our work there 
has been no sign of what one traditionally expects to find in a microscopic description of 
superconductivity: a discussion of a pairing mechanism and a Lagrangian for the degrees 
of freedom that form the Cooper pairs. 

Recall that Landau-Ginzburg theory is the effective field theory description of super- 
conductors near the superconducting phase transition. The dynamical degree of freedom 
is the complex order parameter, <p = {O), that is coupled to a background electromagnetic 
field. The free energy density in the order parameter is 

A/l-G = 2^rl(V + iqA)v\ 2 + aM 2 + ^| 4 • (8.1) 

In this expression m*,q,a and b are phenomenological parameters. The quantity a goes 
through zero at the critical temperature a ~ a(T — T c ). From this free energy one can derive 
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important quantities such as the superconducting coherence length, £, and the superfluid 
density 

* ~ , !u/ 2 ' n « ~ I ■ ( 8 - 2 ) 

Landau-Ginzburg theory is not a microscopic theory. If we wished to study low tem- 
peratures, away from the critical temperature, we would have to supplement the above 
expression with an infinite number of coefficients describing a general functional of ip. The 
curve |</?(T)|, for instance, is an input to rather than an output from this functional. For 
that one needs BCS or some other microscopic theory. The usefulness of Landau-Ginzburg 
theory near the critical temperature is that it relates various experimental quantities and 
can describe the interaction of a superconducting condensate with an electromagnetic field. 

The structural similarity with our phenomenological holographic theory is that we also 
have an infinite number of undetermined parameters at low temperatures. We chose a 
potential with only a mass term, but we could have chosen an arbitrary function of ip. We 
could also have taken a nonminimal coupling between the scalar and the Maxwell field. 
Although this question remains to be fully investigated, we believe that our numerical 
results for (0(T)), for instance, will depend significantly on the gravitational action. As with 
Landau-Ginzburg theory, near the critical temperature our model is much more constrained 
- only the mass term is important to lowest order. 

Despite these similarities, there are three important differences. First, the instability 
which leads to the superconducting phase transition in the CFT has a more natural inter- 
pretation in the gravity theory than in Landau-Ginzburg theory. Gravitationally, given a 
charged black hole and a fixed mass scalar field, the scalar will typically develop a nontrivial 
profile at sufficiently low temperature. Moreover, the curvature of the geometry stabilizes 
the instability without need for higher order terms in the scalar potential. In contrast, in 
Landau-Ginzburg theory, a temperature dependent mass term is added by hand and then 
stabilized by an additional quartic interaction. 

Second, there is a natural way to promote our phenomenological holographic super- 
conductor into a full microscopic description: If we had realised our model as a limit of 
string theory, then the potential for ip would be completely fixed and there would be no 
free parameters. We would have a concrete CFT that underwent a superconducting phase 
transition at a critical temperature specified by the background charge density. Further- 
more, in this theory, the AdS/CFT correspondence allows us to compute all the quantities 
for this superconductor which would normally follow from a BCS-like treatment: the gap 
as a function of temperature, the frequency dependent conductivity, the magnetic pene- 
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tration depth, etc. We have shown how to use AdS/CFT to compute these quantities in 
this paper and we see that the 'feel' of the computation is completely different from weakly 
coupled BCS-like theories. Nonetheless, AdS/CFT applied to a model embedded in string 
theory would be an honest-to-goodness microscopic computation of these quantities in a 
well-defined theory. 

Thirdly, the physical meaning of the potential appearing in the bulk description is com- 
pletely different from the potential in Landau- Ginzburg theory. The AdS/CFT dictionary 
repackages the degrees of freedom of the CFT to make manifest the classicality of a large 
N limit. The theory remains strongly coupled. AdS/CFT is not effective field theory. Our 
phenomenological choice of a 'minimal' model is not guided by Wilsonian arguments but 
rather by simplicity in terms of the degrees of freedom arising through the AdS / CFT dic- 
tionary. This approach makes sense if one accepts the AdS/CFT correspondence as the 
natural tool for an analytic description of strongly coupled theories. The ultimate test of 
this assumption will be the success and robustness of predictions from phenomenological 
AdS/CFT in modeling superconductivity in experimental systems where strong coupling 
and perhaps scale invariance (quantum criticality) play a key role. The heavy fermion 
compounds come to mind as good candidate systems [28 . 



9 Summary 

The main points we have made in this paper are as follows: 

• A minimal AdS/CFT superconductor has a bulk description with a metric, a Maxwell 
field and a charged scalar field (section [2]) . 

• If the CFT is placed at a finite charge density and if the scalar is sufficiently light 
and/or sufficiently charged, then there is a charged condensate in the theory below 
a critical temperature, T < T c . We noted there are two distinct reasons why this 
condensation can happen (section [3J . 

• We computed the charged VEV as a function of temperature (figure [T]) and the critical 
temperature as a function of the charge q of the operator that condenses (figure [2| . 

• We computed the frequency dependent conductivity a(oo) (figures [3] and [4]) and showed 
that a gap opens up for T < T c . There is a delta function in Re [a] at lo = corre- 
sponding to an infinite DC conductivity. We also obtained thermal and thermoelectric 
conductivities (section 4) . 
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• We studied the effect of adding a magnetic field to the holographic superconductor and 
argued that they are always type II. Superconducting droplets form as the magnetic 
field is lowered (section 5). 

• We showed that holographic superconductors generate screening currents obeying the 
London equation. A photon coupled to the superconductor acquires a mass (sections 
6 and 7). 

• Using gravity to study superconductivity has a superficial similarity to a Landau- 
Ginzburg description. Key differences include the fact that the phase transition does 
not have to be put in by hand, and a complete microscopic theory could be obtained 
by realizing a model similar to ours as a limit of string theory (section 8). 

The are many remaining questions to address. Besides technical issues such as un- 
derstanding the zero temperature limit better, three pressing directions of research might 
be emphasized: Firstly to find robust ('universal') results from phenomenological holo- 
graphic superconductors and to understand the extent to which these are useful results for 
experimental systems involving nonconventional superconductivity; secondly, to obtain a 
microscopically understood model by embedding a holographic superconductor into string 
theory; and thirdly to find new mechanisms for superconductivity in the AdS/CFT corre- 
spondence, perhaps with d-wave symmetry or where the critical temperature T c is set by a 
dynamical scale A rather than a charge density. 
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A Instability with a neutral scalar field 

In this appendix we give a proof of the following claim made in the text: When coupled 
to a neutral scalar field with m 2 = —2, the Reissner-Nordstrom AdS black hole becomes 
unstable near extremality. We shall prove this using test functions and the Rayleigh-Ritz 
method. 
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Let us write the black hole metric in terms of the coordinate z = 1/r 

(A.l) 



ds 2 = \ 



dz l 

-f{z)dt 2 + — + dx 2 + dy 2 
J\ z ) 



where 

/ = 1 - (1 + c 2) z 3 + C 2 Z 4 _ (A 2) 

Here c denotes a dimensionless charge density, obtained by rescaling the horizon to z = 1. 
It is related to the physical charge density p and temperature T by [25] 

p_ 16vr 2 c 

T 2 (3 - c 2 ) 2 ' 1 J 

We now want to write the equation of motion for the neutral scalar field tp in Schrodinger 
form. This form will help us gain intuition about the behaviour of the field. The rewriting 
requires rescaling the field and changing variables from z to a new coordinate s. Let 

ds 1 / . . \ 

d~z = r (A - 4) 

(Although irrelevant for the following argument, it follows that the range of s is from 
(the boundary) to oo (the horizon)). Then, taking ^f(s,t) with a time dependence e~ lujt , 
we obtain the Schrodinger equation 

"S" + V(s)* = lu 2 *, (A.5) 
with the potential, written in terms of the z{s) variable, 

y=-f(h + i-%)- <«» 

If this Schrodinger equation has a negative energy bound state, then we have found an 
instability of the black hole. Negative energy in this context implies that uj 2 < 0. Hence u 
is pure imaginary and there are solutions which grow exponentially in time. 



This potential (A. 6) is positive everywhere unless the charge density c is close to the 
extremal value c = v3 (i.e. T = 0). More specifically, the potential develops a negative 
region in the vicinity of the horizon for c > 1. So any instability is restricted to the charge 
values 1 < c < y/3, i.e. 4vr 2 < p/T 2 < oo. 

We cannot solve this Schrodinger equation exactly. We solved it numerically in the main 
text. However, to show the existence of a negative energy bound state it is sufficient to find 
a test function, satisfying the correct boundary conditions, which gives a negative energy. 
The action to use depends on the boundary conditions of the field The general allowed 
falloff at the boundary z — > is 

* ~ a + bz H . (A.7) 
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We can consider either the first or the second of these terms to be the 'non-normalisable' 
mode. The action must be stationary under variations of the normalisable mode. If we 
impose the boundary condition 5a = (i.e. (^(O) = 0) the following action is stationary 
on solutions to the Schrodinger equation 



Ssa=0 



ds 



ds) + ^ 



(A.8) 



However, if we wish to impose 5b = (i.e. 5^'(0) = 0) then we must add a boundary term 

2 



Ssb= 







ds 



d^Y 
~ds~) 



+ {V(s)-oj 2 )^ 



24»- 



d^ 

ds 



(A.9) 



Both of these actions are finite on solutions to the Schrodinger equation, partly due to a 



cancellation in the potential (A. 6) as z — ► which only occurs at the mass we have chosen, 

2 



rn 



-2. In fact, the boundary term in (A.9) vanishes on shell for 'normalisable' modes. 



To show an instability, we need to find test functions such that S(uj = 0) < 0. Let us 
start with the second of the boundary conditions above. Normalisable modes therefore have 
b = 0. A simple test function that satisfies this boundary condition is 



test 



1 



az 



It is easy to check that the (cj = 0) action is minimised by 

21(3c 2 - 5) 



a 



(A.10) 



(A.ll) 



10(9c 2 - 35) ' 

This function then leads to an action which is negative for 1.609 < c < y/3 ~ 1.732. Note 
that in this expression the lower bound in c is an analytic result, it can be expressed in 
terms of ratios of square roots. If we take a more sophisticated test function we can lower 



the bound a little. For instance, taking a fourth order polynomial in z instead of (A. 10) 
leads to the lower bound c ~ 1.584 which is close to actual value found numerically in the 
main text (c ~ 1.582). Therefore the test function method not only indicates the existence 
of unstable black holes for this boundary condition, but also gives a good estimate of the 
minimal p/T 2 at which the normal phase is unstable. 

We can also consider the boundary condition in which 'normalisable' modes have a = 0. 
For this case we have not found a test function indicating the existence of an instability. 
Naively speaking, this is because the falloff ~ z rather than \& ~ 1 forces more kinetic 
energy into the field. Our full numerics in the main text suggest (or at least, are consistent 
with the idea) that there should be an unstable mode in this case also, but that the critical 
value of c should be very close to the extremal value c = a/3- 
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B Analytic estimate of n s for (O2) case 

In this appendix, we attempt to calculate n s analytically at low temperatures for the dimen- 
sion two case, using the method in section (6.1). If we assume that at very low temperature, 



V / 2V' ~ (£ > 2)/r 2 , then (6.5) becomes 

{to 2 -k 2 -q 2 (O 2 } 2 z 2 )A x +A x = . (B.l) 

where z = 1/r and a dot denotes d/dz. This differential equation can be solved in terms of 
parabolic cylinder functions, D v {cz) where the choice 

2 , ,2 



" = -l+ k ^ and c=^aW) (B.2) 

gives the proper exponential fall-off as z gets large. Here the condition that k and u are 
small is more precisely k 2 , to 2 <C q(02). Expanding D u (cz) near the boundary, we find 

A x = ajl- ^gMy^) z + otffj f (b.3) 

where we have suppressed corrections in (k 2 — uj 2 )/q{02)- The London equation here is 
then 

2T(3/4) 



J x = -^^VoW)a x (B.4) 

Numerically, this estimate of n s appears to be wrong by about 25% at low temperatures. 
While 2r(3/4)/r(l/4) w 0.676, the real constant of proportionality appears to be about 
0.546. 
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